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Abstract 

We study the combinatorial FIFO Stack-Up problem, where bins have to be stacked-up 
from conveyor belts onto pallets. This is done by palletizers or palletizing robots. Given k 
sequences of labeled bins and a positive integer p, the goal is to stack-up the bins by iteratively 
removing the first bin of one of the k sequences and put it onto a pallet located at one of p 
stack-up places. Each of these pallets has to contain bins of only one label, bins of different 
labels have to be placed on different pallets. After all bins of one label have been removed 
from the given sequences, the corresponding stack-up place becomes available for a pallet 
of bins of another label. All bins have the same size. The FIFO Stack-Up problem asks 
whether there is some processing of the sequences of bins such that at most p stack-up places 
are used. 

In this paper we strengthen the hardness of the FIFO Stack-Up shown in [20] by consid¬ 
ering practical cases and the distribution of the pallets onto the sequences. We introduce a 
digraph model for this problem, the so called decision graph, which allows us to give a breadth 
first search solution of running time 0{n^ ■ (m-|-2)*^), where m represents the number of pallets 
and n denotes the total number of bins in all sequences. Further we apply methods to solve 
hard problems to the FIFO Stack-Up problem. Therefor we consider restricted versions of 
the problem, two integer programming models, exponential time algorithms, parameterized 
algorithms, and approximation algorithms. In order to evaluate our algorithms, we intro¬ 
duce a method to generate random, but realistic instances for the FIFO Stack-Up problem. 

Our experimental study of running times shows that the breadth first search solution on the 
decision graph combined with a cutting technique can be used to solve practical instances 
on several thousands of bins of the FIFO Stack-Up problem. Further we analyze our two 
integer programming approaches implemented in CPLEX and GLPK. As expected CPLEX 
can solve the instances much faster than GLPK and our pallet solution approach is much 
better than the bin solution approach. 
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1 Introduction 


This paper is about the combinatorial problem of stacking up bins from conveyor belts onto 
pallets. This problem originally appears in stack-up systems that play an important role in delivery 
industry and warehouses. Stack-up systems are often the back end of order-picking systems. 
A detailed description of the applied background of such systems is given in [3 [26]. Logistic 
experiences over 30 years lead to high flexible conveyor-based stack-up systems in delivery industry. 
We do not intend to modify the architecture of existing systems, but try to develop efficient 
algorithms to control them. 

The bins that have to be stacked-up onto pallets reach the stack-up system on a main conveyor 
belt. At the end of the conveyor belt the bins enter the palletizing system. Here the bins are 
picked-up by stacker cranes or robotic arms and moved onto pallets, which are located at stack-up 
places. Often vacuum grippers are used to pick-up the bins. This picking process can be performed 
in different ways depending on the architecture of the palletizing system (single-line or multi-line 
palletizers). Full pallets are carried away by automated guided vehicles, or by another conveyor 
system, while new empty pallets are placed at free stack-up places. 

The developers and producers of robotic palletizers distinguish between single-line and multi- 
line palletizing systems. Each of these systems has its advantages and disadvantages. 

In single-line palletizing systems there is only one conveyor belt from which the bins are picked- 
up. Several robotic arms or stacker cranes are placed around the end of the conveyor. We model 
such systems by a random access storage which is automatically replenished with bins from the 
main conveyor, see Figure |T| The area from which the bins can be picked-up is called the storage 
area. The storage area is the last part of the conveyor belt where the stacker cranes or robotic 
arms reach the bins. 


storage area main conveyor belt 
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stack-up places 


Figure 1: The single-line stack-up system using a random access storage of storage capacity 5. 
The colors represent the pallet labels. Bins with different colors have to be placed on different 
pallets, bins with the same color have to be placed on the same pallet. 


In multi-line palletizing systems there are several buffer conveyors from which the bins are 
picked-up. The robotic arms or stacker cranes are placed at the end of these conveyors. Here, 
the bins from the main conveyor of the order-picking system first have to be distributed to the 
multiple infeed lines to enable parallel processing. Such a distribution can be done by some cyclic 
storage conveyor, see Figure |2l From the cyclic storage conveyor the bins are pushed out to the 
buffer conveyors. A stack-up system using a cyclic storage conveyor is, for example, located at 
Bertelsmann Distribution GmbH in Giitersloh, Germany. On certain days, several thousands of 
bins are stacked-up using a cyclic storage conveyor with a capacity of approximately 60 bins and 
24 stack-up places, while up to 32 bins are destined for a pallet. This palletizing system has 
originally initiated our research. 

If we ignore the task to distribute the bins from the main conveyor to the k buffer conveyors, 
i.e., if the buffer conveyors are already filled, and if each arm can only pick-up the first bin of one 
of the buffer conveyors, as it is the case if stacker cranes are used, then the system is called a 
FIFO palletizing system. Such systems can be modeled by several simple queues, see Figure O 

From a theoretical point of view, an instance of the FIFO Stack-Up problem consists of 
k sequences qi,... ,qk of bins and a number of available stack-up places p. Each bin of q is 
destined for exactly one pallet. The EIFO Stack-Up problem is to decide whether one can 
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Figure 2: A multi-line stack-up system with Figure 3: The FIFO stack-up system ana- 

a pre-placed cyclic storage conveyor. lyzed in this paper. 


remove iteratively the bins from the sequences such that in each step only one of the first bins of 
qi,... ,qk is removed and after each step at most p pallets are open. A pallet t is called open, if 
at least one bin for pallet t has already been removed from the sequences, and if at least one bin 
for pallet t is still contained in the remaining sequences. If a bin b is removed from a sequence 
then all bins located behind b are moved-up one position to the front (cf. Section [2] for the formal 
definition). 

Every processing should absolutely avoid blocking situations. A system is blocked, if all stack- 
up places are occupied by pallets, and none of the bins that may be used in the next step are 
destined for an open pallet. To unblock the system, bins have to be picked-up manually and moved 
to pallets by human workers. Such a blocking situation is sketched in Figure [3] 

The single-line stack-up problem can be defined in the same way. An instance for the single- 
line stack-up problem consists of one sequence q of bins, a storage capacity s, and a number of 
available stack-up places p. In each step one of the first s bins of q can be removed. Everything 
else is defined analogously. 

Many facts are known about single-line stack-up systems [551 HZl HH]- In [IS] it is shown 
that the single-line stack-up decision problem is NP-complete, but can be solved efficiently if the 
storage capacity s or the number of available stack-up places p is a fixed constant. The problem 
remains NP-complete as shown in [27] . even if the sequence contains at most 9 bins per pallet. 
In [27j , a polynomial-time off-line approximation algorithm for minimizing the storage capacity s 
is introduced. This algorithm yields a solution that is optimal up to a factor bounded by log(p). 
In [28j the performances of simple on-line stack-up algorithms are compared with optimal off-line 
solutions by a competitive analysis 13 [g. 

The EIFO Stack-Up problem has not been studied by other authors up to now, although 
stack-up systems play an important role in delivery industry and warehouses. In our studies, we 
neither limit the number of bins for a pallet nor restrict the number of stack-up places to the 
number of buffer conveyors. That is, the number of stack-up places can be larger than or less 
than the number of buffer conveyors. We could show in [20] that the FIEO Stack-Up problem 
is NP-complete, but can be solved in polynomial time if the number k of sequences or the number 
p of stack-up places is fixed. 

This paper is organized as follows. In Section |2l we give preliminaries and the problem state¬ 
ment. In Section 13 we recall the definition of a sequence graph from |20] . This digraph has a 
vertex for every pallet and an arc from pallet a to pallet b, if and only if in any processing pallet 
b can only be closed if pallet a has already been opened. We give an algorithm which allows us 
to compute the sequence graph in time 0{n k ■ rn^), where m represents the number of pallets 
and n denotes the total number of bins in all sequences. Further we show that the hardness of 
the FIFO Stack-Up problem even holds for the practical case k < m and we the discuss the 
influence of the distribution of the pallets onto the sequences. 

In Section |3 we consider two further digraph models that allow us to find dynamic program- 
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ming solutions for the FIFO Stack-Up problem. The first digraph is the processing graph. It 
has a vertex for every possible configuration of the system and an arc from configuration A to 
configuration B if configuration B can be obtained from configuration ^ by a single processing 
step. The algorithmic use of the processing graph was already mentioned in |20| and will be 
explained more in detail here in order to ease the understanding on the second digraph model, 
which is called the decision graph. It has a vertex for every decision configuration of the system, 
i.e. for configurations such that for every sequence the next bin is destined for a non-open pallet 
and an arc from decision configuration A to decision configuration B if configuration B can be 
obtained from configuration A by a processing step including automatic steps. The decision graph 
allows us to give a breadth first search solution for the FIFO Stack-Up problem of running time 
• (m + 2)^). 

In Section[5]we apply methods to solve hard problems to the FIFO Stack-Up problem. There¬ 
for we consider restricted versions of the problem, exponential time algorithms, and approximation 
algorithms. We also give two integer programming models to solve the FIFO Stack-Up problem. 
The first model computes a bin solution, i.e. an order in which the bins can be removed, and the 
second one computes a pallet solution, i.e. an order in which the pallets can be opened. Both 
models are using a polynomial number of variables and a polynomial number of constraints to 
compute the minimum number of stack-up places. Further we study the fixed-parameter tractabil- 
ity of the problem. The idea behind fixed-parameter tractability is to split the complexity into 
two parts - one part that depends purely on the size of the input, and one part that depends on 
some parameter of the problem that tends to be small in practice. Based on our three digraph 
models and our two integer programming models we give parameterized algorithms for various 
parameters which imply efficient solutions for the FIFO Stack-Up problem restricted to small 
parameter values. 

In Section ini we introduce a method to generate random, but realistic instances for the FIFO 
Stack-Up problem. We generated instances on several thousand bins of the FIFO Stack-Up 
problem which could be solved by our breadth first search solution combined with some cutting 
technique on the decision graph in a few minutes. Further we analyze two integer programming 
approaches implemented in CPLEX and GLPK. As expected CPLEX can solve the instances much 
faster than GLPK and our pallet solution approach is much better than the bin solution approach. 


2 Problem Statement 

We consider sequences 

Ql — (^1; • ■ • ; . . . , + ■ • ■ ; ■ ■ ■ : Qk i 

of bins. All these bins are pairwise distinct. These sequences represent the buffer queues (handled 
by the buffer conveyors) in real stack-up systems. Each bin b is labeled with a pallet symbol plt(b) 
which can be some positive integer. We say bin b is destined for pallet pltib). The labels of the 
pallets can be chosen arbitrarily, because we only need to know whether two bins are destined for 
the same pallet or for different pallets. The set of all pallets of the bins in some sequence qi is 
denoted by 

plts{qi) = {plt{h) I b e qi}. 

For a list of sequences Q = (qi,... ,qk) we denote 

plts{Q) = plts{qi) U • • • U plts[qk). 

For some sequence g = (6i,..., 6„) we say bin bi is on the left o/bin bj in sequence g if i < j. And 
we say that such a bin bi is on the position i in sequence g, i.e. there are i — 1 bins on the left 
of b in sequence g. The position of the first bin in some sequence qi destined for some pallet t is 
denoted by first{qi,t), similarly the position of the last bin for pallet t in sequence qi is denoted 
by last{qi,t)- For technical reasons, if there is no bin for pallet t contained in sequence qi, then 
we define first{qi,t) = \qi\ 1, and last{qi,t) = 0. 
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Let Q = (< 71 ,..., gfe) be a list of sequences, and let Cq = {ii,... ,ik) be some tuple in Such 
a tuple Cq is called a configuration, if 0 < ij < \qj\ for each sequence qj € Value ij denotes 
the number of bins that have been removed from sequence qj, see Example 12.21 

A pallet t is called open in configuration Cq = {ii,..., ik), if there is a bin for pallet t at some 
position less than or equal to ij in sequence qj, i.e. first{qj,t) < ij, and if there is another bin 
for pallet t at some position greater than i^ in sequence qe, i.e. last(qi,t) > ii, see Example 12.21 
In view of the practical background we only consider sequences that contain at least two bins for 
each pallet. The set of open pallets in configuration Cq is denoted by open{CQ), the number of 
open pallets is denoted by ffopeniJjQ). 

Remark 2.1 Within several algorithms we will need the set of open pallets within some config¬ 
uration Cq. Therefore we first compute all the values first{qi,f) and last{qi,t) for qi € Q and 
t G plts{Q) in time 0{k ■ max{|gi |,... ,\qk\} k ■ m) respectively 0(n + fc ■ m). Using these values 
we ean test in time 0{k) whether some pallet t is open within a given eonfiguration. By performing 
this test for eaeh of the m pallets, for some eonfiguration Cq, we ean compute ffopen{CQ) in time 
Oim-k). 

A pallet t G plts{Q) is called elosed in configuration Cq, if last{qj,t) < ij for each sequence 
qj G Q. Initially all pallets are unproeessed. From the time when the first bin of a pallet t has 
been removed from a sequence, pallet t is either open or closed. 

For some configuration (ii,... ,i}f) we define 

front((ii,..., ik)) = {plt{b) \ I < j < fc, 6 is on position ij + I in sequence qj}. 

Informally speaking front((ii,... ,ik)) is the set of all pallets of the first bins of the remaining 
sequences in configuration (ii,..., ik). (cf. Example 12.21 and Table[T|). 

Let Cq = {ii,... ,ik) be a configuration. The removal of the bin on position ij + I from 
sequence qj is called a transformation step. A sequence of transformation steps that transforms 
the list Q oi k sequences from the initial configuration (0,0,... ,0) into the final configuration 
(kill | 92 |j ■ • ■ j |<?fe|) is called a processing of Q, see Example 12.21 

It is often convenient to use pallet identifications instead of bin identifications to represent a 
sequence q. For r not necessarily distinct pallets t\,... ,tr let [ti,..., tr] denote some sequence of 
r pairwise distinct bins (hi,..., hr), such that plt(bi) = ti iov i = 1,... ,r. We use this notation for 
lists of sequences as well. Furthermore, for some positive integer value n, let [n] := {1,2,..., n} 
be the set of all positive integers between 1 and no 

Example 2.2 (Processing) Consider the list Q = (( 71 ,( 72 ) of two sequences 

qi = (hi,...,hi) = [a,b,a,b] 

and 

(72 = (& 5 , • ■ •, & 10 ) = [c, d, c, d, a, 6]. 

TableUi shows a proeessing of Q with 2 staek-up places. The underlined bin is always the bin that 
will be removed in the next transformation step. The already removed bins are shown greyed out. 

We consider the following problem. 

Name: FIFO Stack-Up 

Instance: A list Q = (qi,... ,qk) oi k sequences of bins, for every bin b oi Q its pallet symbol 
plt(b), and a positive integer p. 

Question: Is there a processing of Q, such that in each configuration during the processing of Q 
at most p pallets are open? 

^An alternative definition of configurations using subsequences was given in ED]. 

^We will use square brackets in several different notations. Although the meaning becomes clear from the context 
we want to emphasize this fact. 


5 



i 

qi 

92 

Cq 

front {Cq) 

open{CQ) 

bin to remove 

0 

[a, b, a, b] 

[c, d, c, d, a, b] 

( 0 , 0 ) 

{a,c} 

0 

bb 

1 

[a, b, a, b] 

[c.d, c, d, a, b] 

( 0 , 1 ) 

{a,d} 

{4 

be 

2 

[a, b, a, b] 

[c,d.c, d,a,b] 

( 0 , 2 ) 

{a,c} 

{c,d} 

67 

3 

[a, b, a, b] 

[c, d, c. d, a, b] 

(0,3) 

{a,d} 

{d} 

bs 

4 

[a, b, a, b] 

[c, d, c, d, a, b] 

(0,4) 

{a} 

0 

bi 

5 

[a. b, a, b] 

[c, d, c, d, a, b] 

(1,4) 

{a,b} 

{a} 

bg 

6 

[a,b, a, b] 

[c, d, c, d, a. 6 ] 

(1,5) 

{b} 

{a} 

b2 

6 

[a, b, a, b] 

[c,d,c, d,a,b] 

(2,5) 

{a,b} 

{a,b} 

bio 

7 

[a, b. a, b] 

[c, d, c, d, a, b] 

( 2 , 6 ) 

{a} 

{a,b} 

bs 

8 

[a, b, a, 5] 

[c, d, c, d, a, b] 

(3,6) 

{b} 

{b} 

bi 

9 

[a, b, a, b] 

[c, d, c, d, a, b] 

(4,6) 

0 

0 

- 


Table 1: A processing of Q = (<Zij 92 ) from Example 12.21 with 2 stack-up places. In this simple 
example it is easy to see that there is no processing of Q that needs less than 2 stack-up places, 
because pallets a and b as well as c and d are interlaced. 


We use the following variables in the analysis of our algorithms: k denotes the number of 
sequences, and p stands for the number of stack-up places, while m represents the number of pallets 
in plts{Q), and n denotes the total number of bins, i.e. n = Uk- Finally, N = max{|( 7 i|,..., |gfe|} 
is the maximum sequence length. 

For some instance / of the FIFO Stack-Up problem numbers are encoded binary and se¬ 
quences are encoded by sequences of pallet symbols thus the size |/| can be bounded by 

|/| e C>(n-log2(m)-f log2(p)). (1) 

The size of the input is important for the analysis of running times in Section [5l 

Remark 2.3 In view of the practical background, it holds p < m, otherwise each pallet eould be 
staeked-up onto a different stack-up place. Furthermore, k < m, otherwise all bins of one pallet 
could be channeled into one buffer queue in the multi-line stack-up systems with pre-placed cyclic 
storage conveyor, see Figure [3 Finally m < ^ < n, since there are at least two bins for each 
pallet. 

The relation n < k ■ N and the assumption m < ^ imply the following bound. 

Corollary 2.4 m < i.e. m G 0{k ■ N), 

The following estimation can be shown by induction on k. 

Corollary 2.5 k ■ N < {N -\- 1)^ 

For the solution of the FIFO Stack-Up problem for some list of sequences Q we use the 
notation of a bin solution and of a pallet solution, which are defined as follows. Let B = 
(& 7 r(i)j ■ • ■ j & 7 r(n)) be the order in which the bins are removed during a processing of Q. Then 
B is called a bin solution of Q. In Example 12.21 we have 

B = {b5,bG,b7,bs,bi,bQ,b2,bio,b3,b4) = [c,d,c,d,a,a,b,b,a,b] 


as a bin solution. 

Let T = (ti,..., tm) be the order in which the pallets are opened during the processing of Q. 
Then T is called a pallet solution of Q. In Example 12.21 we have 

T = (c, d, a, b) 


as a pallet solution. 
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3 NP-hardness 


Next we recall the connection between the used number of stack-up places for a processing of an 
instance Q and the directed pathwidth of the sequence graph Gq from [50] , which is useful for our 
hardness results, our integer programming model for computing a pallet solution, and for several 
parameterized algorithms. 

3.1 Directed Pathwidth 

According to Barat [T], the notion of directed pathwidth was introduced by Reed, Seymour, 
and Thomas around 1995 and relates to directed treewidth introduced by Johnson, Robertson, 
Seymour, and Thomas in [21]. A directed path-decomposition of a digraph G = {V,A) is a 
sequence (Xi,..., Ar) of subsets of V, called bags, such that the following three conditions hold 
true. 

(1) XiU...UXr = V, 

(2) for each (u,v) G A there is a pair i < j such that u G Xi and v G Xj, and 

(3) for all i,j, £ with l<i<j<£<r it holds XifiXe C Xj. 

The width of a directed path-decomposition X = (Xi,... ,Xr) is maxi<i<r |Xi| — 1. The directed 
pathwidth of G, d-pw(G) for short, is the smallest integer w such that there is a directed path-de¬ 
composition for G of width w. For symmetric digraphs, the directed pathwidth is equivalent to the 
undirected pathwidth of the corresponding undirected graph |24| , which implies that determining 
whether the pathwidth of some given digraph is at most some given value w is NP-complete. For 
each constant w, it is decidable in polynomial time whether a given digraph has directed pathwidth 
at most w, see Tamaki [29]. 

3.2 The Sequence Graph 

The sequence graph Gq = {V,A) for an instance Q = {qi, ■ ■ ■ ,qk) is defined by vertex set V = 
plts{Q) and the following set of arcs. There is an arc {u,v) G A if and only if there is a sequence 
qi = (bni_i+i, ■. •, bni) with two bins bi, bj such that i < j, plt(bi) = u, plt(bj) = v, and u ^ v. 

Example 3.1 (Sequence Graph) Figure\^shows the sequence graph Gq for Q = ( 91 , 92 , 93 ) 
with sequences 91 = [a, a, d, e, d], 92 = [c, b, b, d], and 93 = [c, c, d, e, d]. 



Figure 4: Sequence graph Gq of Example 13. II 

If Gq = {V,A) has an arc {u,v) G A then u ^ v and for every processing of Q, pallet u is 
opened before pallet v is closed. Digraph Gq = {V, A) can be computed in time 0{n + k ■ mf) by 
the algorithm Create Sequence Graph shown in Figure [5] 

A value is added to a list only if it is not already contained. To check this efficiently in time 
0(1) we have to use a boolean array. In our algorithm V and L are implemented by boolean 
arrays. Therefore we need some preprocessing phase where we run through each sequence and 
seek for the pallets. This can be done in time 0{n + k ■ m) C 0{n + m^). 

In [20] we have shown the following correlation between the used number of stack-up places 
for a processing of an instance Q and the directed pathwidth of the sequence graph Gq . 
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Algorithm Create Sequence Graph 


for each sequence q G Q do 

b := first bin of sequence q 

add plt{b) to vertex set V, if it is not already contained 
L := {plt{b)) l> L contains pallets of bins up to bin b 

for i := 2 to Ig] do 

b := i-th bin of sequence q 

add plt{b) to vertex set V^ if it is not already contained 
if {i = last{q, pU{b))) 

for each pallet t G L do 

a t ^ plt(b) add arc (t, plt(b)) to arc set A, if it is not already contained 
if {i = first{q,plt{b))) 
append{plt{b), L) 


Figure 5: Create the sequence graph G = iV.A) for some given list of sequences Q. 

Theorem 3.2 Let Q = (gi,..., qk) then digraph Gq = {V, A) has directed pathwidth at most p— 1 
if and only if Q can be processed with at most p stack-up places. 

This characterization extends the previously known areas of applications for directed pathwidth 
in graph databases and boolean networks, which have been shown in [^. 

3.3 Hardness Result 

Next we will discuss the hardness of the FIFO Stack-Up problem. In contrast to Section !??^ we 
transform an instance of a graph problem into an instance of the FIFO Stack-Up problem. 

Let G = {V,A) be some digraph and A = {ai,...,af} its arc set. The sequence system 
Qc = {qi,. ■ ■ ,qe) for G is defined as follows. 

(1) There are 2£ bins &i,..., 62 ^. 

(2) Sequence qt = (& 2 i-i, ^ 2 i) for 1 < * < .^. 

(3) The pallet symbol of bin 621-1 is the first vertex of arc Oi and the pallet symbol of 621 is the 
second vertex of arc Oi for I < i < i. Thus plts{QG) = V. 

Example 3.3 (Sequence System) For the digraph G of Figure the corresponding sequence 
system is Qg = (gi, 92 , 93 , < 74 , <? 5 , 96 , 9 ?); where 

9 i = [a,&], 92 = [ 6 ,c], ga = [c,d], q^ = [d,e], q^ = [e,a], 

96 = [e,/], 97 = [f,a]. 

The sequence graph of Qg is G. 



Figure 6: Digraph G of Example 13.31 

The definition of sequence system Qg and sequence graph Gq, defined in Section [3^ imply 
the following results shown in |20) . 







Proposition 3.4 ( |20p For every digraph G it holds G = Gq^. 

Theorem 3.5 ([2Q1) The FIFO Stack-Up problem is NP-complete, even if the sequences of Q 
contain together at most 6 bins per pallet. 

Next we strengthen the hardness by considering further restricted versions of the FIFO Stack- 
Up problem. First we can bound the maximum sequence length N, since by the definition of the 
sequence system we obtain instances where iV = 2 in the proof of Theorem 13.51 

Corollary 3.6 The FIFO Stack-Up problem is NP-complete, even if N is bounded by some 
constant greater than 1. For N = 1 the FIFO Stack-Up problem can obviously be solved in 
polynomial time. 

In order to consider the distribution of the bins of a pallet t onto the sequences we define 

dgit) = \{q£Q \ t£ plts{q)}\ 


and 


dq = maxdQ(t). 

tG[m] 


By Theorem 13.51 we have shown the following result. 


Corollary 3.7 The FIFO Stack-Up problem is NP-complete, even if dg = 6. For dg = 1 the 
FIFO Stack-Up problem can be solved in polynomial time, since we can process all sequences 
one after the other. 

For dg € {2,..., 5} the complexity of the FIFO Stack-Up problem remains open. 

In Remark l2.3l we have restricted to practical instances where k < m. This is not given within 
the hardness results of |20) . 

Corollary 3.8 The FIFO Stack-Up problem is NP-complete, even if k < m and the sequences 
of Q contain together at most 6 bins per pallet. 


Proof In order to carry over the hardness to k < m, we can modify a given list of sequences 
Q = (gi,..., qk) as follows. We introduce 3fc new pallets at, bi, a ior 1 < i < k. We define new 
sequences qi, 1 < i < k, from the old ones by 

qi — qi [oi j o,ijbi,bi, Ci, Ci], 

i.e. by concatenating qi and [ai,ai,bi,bi,Ci,Ci]. Since list Q can be processed with at most p 
stack-up places, if and only if list Q' = {q[,... ,q'^f) can be processed with at most p stack-up 
places and the number of pallets in Q' can be increased arbitrarily, the FIFO Stack-Up problem 
is also hard for the case k < m. □ 


□ 


4 Dynamic Programming Algorithms to Solve the FIFO 
Stack-Up Problem 

Our aim in controlling FIFO stack-up systems is to compute a processing of the given sequences of 
bins with a minimum number of stack-up places. Such an optimal processing can always be found 
by computing the processing graph or the decision graph. The algorithmic use of the processing 
graph was already mentioned in |20) and will next be explained in more detail in order to ease the 
understanding in the subsequent section on the decision graph. 
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4.1 The Processing Graph 

The processing graph G = {V, A) contains a vertex for every possible configuration. Each vertex 
V representing some configuration Cq[v) is labeled by the number ^open^Cq^v)). There is an 
arc from vertex u representing configuration (ui,... ,Uk) to vertex v representing configuration 
(til ,... ,Vk) if and only if = Vi — 1 for exactly one element of the vector and for all other elements 
of the vector uj = Vj. The arc is labeled with the bin that will be removed in the corresponding 
transformation step. 

Example 4.1 (Processing Graph) For the sequences of Example \2.S\ we get the processing 
graph of Figure 



Figure 7: The processing graph of Example 12.21 Instead of the bin each arc is labeled with the 
pallet symbol of the bin that will be removed in the corresponding transformation step. The 
shaded vertices will be important in the following section. The upper left vertex represents the 
initial configuration and the lower right vertex represents the final configuration. 

Obviously, every processing graph is directed and acyclic. Every bin solution describes a 
path from the initial configuration ( 0 , 0 ,... , 0 ) to the final configuration (|gi|, \q 2 \^.. ■, \qk\) in the 
processing graph. We are interested in such paths where the maximal vertex label on that path 
is minimal. 

The processing graph can be computed in time 0{k ■ {N + 1)^) by some breadth first search 
algorithm as follows. We store the already discovered configurations, i.e. the vertices of the graph, 
in some list L. Initially, list L contains only the initial configuration. In each step of the algorithm 
we take the first configuration out of list L. Let Cq = {ii,... ,ik) be such a configuration. Eor 
each j G [fc] we remove the bin on position ij + I from sequence qj , and get another configuration 
Cq. We append Cq to list L, add it to vertex set E, if it is not already contained, and add 
(CqjC'q) to arc set A. 

Eor each configuration we want to store the number of open pallets of that configuration. This 
can be done efficiently in the following way. Eirst, since none of the bins has been removed from 
any sequence in the initial configuration, we have ffopen{(f), . .., 0)) = 0. In each transformation 
step we remove exactly one bin for some pallet t from some sequence q^, thus 

#ope?r((fi,...,fj-i,fj+ l,fj+i,...,4)) = #open((fi,...,fj_i,'ij,'ij+i,...,ffe)) + Cj (2) 

where Cj = I if pallet t has been opened in the transformation step, and Cj = —I if pallet t has 
been closed in the transformation step. Otherwise, Cj is zero. If we put this into a formula we get 

{ I, if first{qj,t) = + 1 and first{q£, t) > ii ^ j 

— 1 , if last{qj,t) = + 1 and last{q£, t) <ii ^ f 

0 , otherwise. 
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Please remember our technical definition of first{q,t) and last{q,t) from page S] for the case that 
t ^ plts{q). 

That means, the calculation of value ^open{CQ{v)) for the vertex v representing configuration 
Cq(v) can be done in time 0{k) if the values first{qj, t) and last{qj, t) have already been calculated 
in some preprocessing phase. Such a preprocessing can be done in time 0{k ■ N + k ■ m) due to 
Remark l2.11 which can be bounded by 0{k ■ (N + 1)^') due to Corollary [23] and Corollarv l2.5l The 
number of vertices is in 0((7V +1)^'), so the vertex labels can be computed in time 0{k ■ {N + 1)^). 
Since at most k arcs leave each vertex, the number of arcs is in 0{k ■ {N + 1)^), and each arc can 
be computed in time 0{1). Thus, the computing can be done in total time 0{k ■ {N + 1)^). 

Let s be the vertex representing the initial configuration, and let topol : V —>■ N be a topological 
ordering of the vertices of the processing graph G = (P, A) such that topol(u) < topol{v) holds for 
each (tt, v) G A. For some vertex v G V and some path P = (vi,..., vg) with vi = s, vg = v and 
G A we define 

valp(v) := open(C q{u))) 

uGP 

to be the maximum vertex label on that path. Let Viv) denote the set of all paths from vertex s 
to vertex v. Then we define 

vallv) := min (valpiv)). 

P&V{v) 

The problem is to compute the value val {t) where t is the vertex representing the final configura¬ 
tion. It holds 


val{v) = max{^open {Cq{v)), min (val{u))}, (3) 

{u,v)£A 

because each path P G V{v) must go through some vertex u with (it, v) G A. So we may iterate over 
the preceding vertices of v instead of iterating over all paths. If ^open{CQ{u)) < #open{CQ{v)) 
for all preceding configurations then a pallet must have been opened in the last step to reach 
configuration Cq(v). 


Algorithm Find Path 


o 

II 

l> Computation according to Equation (|3|) 

for each vertex v ^ s in order of topol do 
val[v\ := oo 

for every (it, u) G A do 
if (yal[u] < val[v]) 
val[v] := val[u\ 
path[v] := It 

[> Compute min(„ „)g^(iia^(it)) =: val{v) 

if {val[v] < ^open[v]) 
val[v] := ^open[v] 

t> Compute max{^open(CQ(ii)),i!a/(i;)} =: val{v) 


Figure 8: Finding an optimal processing by dynamic programming. 

The value val{vg) can be computed by Algorithm Find Path given in Figure H] The corre¬ 
sponding path 

P ={vi,...,vg) (4) 

is obtained by path[vg] = vg-i, path[vg-i] = vg- 2 , ■ ■., path[v 2 ] = vg. For the running time we 
observe that a topological ordering of the vertices of digraph G can be found by a depth first search 
algorithm in time 0(|F| + |A|). The remaining work of algorithm Find Path can also be done in 
time 0(|F| + |A|). In the processing graph we have |F| G 0{{N + 1)^), and |A| G 0{k ■ {N + I)^). 

It is not necessary to explicitly build the processing graph to compute an optimal processing, 
we have done it just for the sake of clarity and to enhance understanding. We combine the 


II 





construction of the processing graph with the topological sorting and the path finding by some 
breadth first search algorithm Optimal Bin Solution shown in Figure [51 Algorithm Optimal 
Bin Solution uses the following two operations. 

• head{L) yields the first entry of list L and removes it from L. 

• append{e, L) adds element e to the list L, if e is not already contained in L. 


Algorithm Optimal Bin Solution 


#open[ ( 0 ,..., 0 ) ] :=0 
val[ (0,...,0) ] := 0 
L:=( (0,...,0) ) 
pred[ ( 0 ,..., 0 ) ] := 0 
while L is not empty do 
C := head{L) 
if {pred[C] ^ 0) 

Extend Path(C') 
for j := 1 to A: do 

Cs . {i \ ^ ^ i j 

if {Cs is not in L) 

compute 41^open[Cs\ according to 
append{Cs, L) 
append(C, pred[Cs\) 


O List of uninvestigated configurations 
O List of predecessors of some configuration 

> let C = (ii,... , 4 ) 

l> all predecessors of C are computed 


Equation ([5]) 


Figure 9: Construction of the processing graph and computation of an optimal bin solution at 
once by breadth first search. 


Algorithm Extend Path(C) 


val[C] := oo 

for each Cp in list pred\C] do 
if {val[C] > val[Cj)]) 
val[C] := val[Cp\ 
path[C] := Cp 
if {val[C] < #opert[C]) 
val[C] := ^open[C] 


\> Compute val[C] due to Equation (|31) 


Figure 10: Submethod to extend a path by one vertex C. 

At the end of the processing of algorithm Optimal Bin Solution the variable path contains 
a path P as shown in (|4]) where the maximal vertex label is minimal. 

The computation of all at most {N +1)^ values ^open{CQ{v)) can be performed in time 
0{k ■ {N + 1)^). A value is added to some list only if it is not already contained. To check this 
efficiently in time 0(1) we use a boolean array over all possible configurations. This array can be 
initialized in time 0{{N + 1)^). Thus, we can compute the minimal number of stack-up places 
necessary to process the given FIFO Stack-Up problem as well as such a bin solution in time 
0{k-{N+ 1)>^). 

Theorem 4.2 The FIFO Stack-Up problem can be solved in time 0{k ■ {N -|- 1)^). 
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4.2 The Decision Graph 

During a processing of a list Q of sequences there are often configurations for which it is easy 
to Hnd a bin b that can be removed such that a further processing with p stack-up places is still 
possible. This is the case, if bin b is destined for an already open pallet. A configuration {ii,... ,ik) 
is called a decision configuration, if the bin on position ij -|- 1 of sequence qj for each j G [fc] is 
destined for a non-open pallet, i.e. front{{ii ,..., ik)) H open((ii ,..., ik)) = 0. We can restrict 
FIFO stack-up algorithms to deal with such decision configurations, in all other configurations the 
algorithms automatically remove a bin for some already open pallet. 

A solution to the FIFO Stack-Up problem can always be found by computing the decision 
graph for the given instance of the problem. The decision graph G = (V, A) has a vertex for 
each decision configuration into which the initial configuration can be transformed. There is 
an arc (u,u) G A from a vertex u representing decision configuration (ui,... ,Uk) to a vertex v 
representing decision configuration (vi,... ,Vk) if and only if there is a bin b on position Uj -|- 1 in 
sequence qj such that the removal of b in the next transformation step and the execution of only 
automatic transformation steps afterwards lead to decision configuration (ui,..., Vk)- Arc (u, v) is 
labeled with the pallet symbol of the bin that will be removed in the corresponding transformation 
step. 

Example 4.3 (Decision Graph) In Figure[T^the decision graph for Examvle \2.2\ is shown. 



Figure II: The decision graph of ExamDle l2.21 It consists of the shaded vertices from the processing 
graph in Figure [T] 

Every decision graph is directed and acyclic. Each pallet solution leads to a path from the 
initial configuration ( 0 , 0 ,... , 0 ) to the final configuration (|gi|, \q 2 \, ■ ■ ■, |gfc|) in the decision graph. 
We are interested in such paths where the maximal vertex label on that path is minimal. 

The decision graph can be computed by some breadth first search algorithm as follows. We 
store the already discovered decision configurations, i.e. the vertices of the graph, in some list 
L. Initially, list L contains only the initial configuration. In each step of the algorithm we take 
the first configuration out of list L. Let Cq = (fi, ... ,ik) be such a decision configuration. For 
each j G [fc] we remove the bin on position ij -|- 1 from sequence qj and execute all automatic 
transformation steps. By this way we reach another decision configuration Cq = {i'l,... ,i'jf) or 
the final configuration, and we put this configuration into list L, we add Cq to vertex set V, and 
we add {Cq,Cq) to arc set A, if it is not already contained. 

We can combine the construction of the decision graph and the path finding by some breadth 
first search algorithm Find Optimal Pallet Solution, as shown in Figure [121 

The running time of algorithm Find Optimal Pallet Solution can be estimated as follows. 
In each sequence qj there are bins for at most m pallets. Each pallet can be opened at most once, 
so in a decision configuration (*i, ... ,ik) it must hold ij -\-l = first[qj,f) for some pallet t. And 
since ij = 0 and ij = \qj\ are possible as well, the decision graph has at most (m -I- 2)^ vertices 
each representing a decision configuration in list L. 

We use a boolean array of size (m -|- 2)^ to indicate whether some decision configuration is 
already contained in L or not. The initialization can be performed in 0{{m + 2)^'). Since there 
are 0{{m + 2 )^') decision configurations each having at most k predecessors we need at most 
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Algorithm Find Optimal Pallet Solution 


open[ (0 ,... ,0) ] := 0 


#open[ (0 ,...,0) ] :=0 


val[ (0,...,0) ] := 0 


L:=( (0,...,0) ) 

> List of configurations to expand 

pred[ (0 ,..., 0) ] := 0 

> List of predecessors of some decision configuration 

while L is not empty do 


C := head{L) 

> let C = (ii,... ,4) 

if {pred[C] ^ 0) 

> all predecessors of C are computed 


Extend Path(C') 
for j -.= 1 to k do 

Cg . “t" 

O := open[C\ iJ pit [hi-j^i) 

for £ := 1 to fc do > perform automatic transformation steps 

while the first bin in qi is destinated for an open pallet in O 

let Ca be the configuration obtained by removing the first bin of qt 
if {Cs is not in L) 

compute 4i^open[Ca\ and open[Cs] according to Remark 12.II 
append{Cs, L) 
append(C, pred[Cs \) 


Figure 12: Construction of the decision graph and computation of an optimal pallet solution at 
once by breadth first search. 


0{k ■ (m + 2)^) tests whether some decision configuration is already contained in L. Since each 
test can be done in time 0{k), we can bound the running time for all tests by 0{k‘^ ■ (to + 2)^) C 
0{n'^ ■ (to + 2 )^). 

The automatic transformation steps can be performed for every j S [k] in time 0{n + k) and 
thus for all decision configurations in L in time 0{{m + 2)^ ■ k ■ {n + k)). The computation of 
4j^open{Ca) for one decision configuration Cs can be done in time 0{m ■ k) by Remark 12.II Since 
we only compute ^open{Cs) for decision configuration Cs which is not already contained in L, for 
all decision configurations we need time in 0{{m + 2 )^ • to • k)). 

Thus the running time of algorithm Find Optimal Pallet Solution is in 0{{m + 2)^ ■ k ■ 
(n + fc) + (to + 2)^ ■ m ■ k)) C 0{{m + 2)^ ■ n^). 

If we have a pallet solution T = (ti,. .. computed by any FIFO stack-up algorithm, we 
can convert the pallet solution into a sequence of transformation steps, i.e. a processing of Q, in 
time 0{n ■ k) C 0{'n?) by some simple algorithm: Repeatedly, in the f-th decision configuration 
choose the pallet ti, remove a bin for pallet U, and execute all automatic transformation steps, 
until the next decision configuration is reached. If no bin for pallet ti can be removed in the f-th 
decision configuration, or if more than p pallets are opened, reject the input. Each transformation 
step of these at most n steps can be done in time 0{k) by Equation ([2)). 

Theorem 4.4 The FIEO Stack-Up problem can be solved in time 0{ii? ■ (to -|- 2)^). 

Another way to describe a path from the initial configuration to the final configuration in the 
decision graph is as follows. Let Si € [fc] denote the sequence from which a bin for pallet ti will be 
removed in the next transformation step after reading the i-th decision configuration on a path. 
Then S = (si,..., Sm) is called a sequenee solution. In Example 12.21 we have S = (2, 2, 1,1) as a 
sequence solution. 

There are at most fc™ sequence orders (si,... ,Sm), where Si G [fc]. The following algorithm 
checks in time 0(n ■ k) C Oin^), whether one sequence order describes a sequence solution: In the 
f-th decision configuration remove the bin from sequence Si. Reject the input, if this is impossible. 
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Otherwise do all automatic transformation steps until the next decision configuration is reached. 
Each step can be done in time 0{k) by Equation (I2|). Reject the input, if more than p pallets are 
opened. 

If we enumerate all possible sequence orders and check each order, we can solve the FIFO 
Stack-Up problem. 

Theorem 4.5 The FIFO Stack-Up problem can be solved in time 0{n^ ■ k™"). 

5 Application of Methods to Solve Hard Problems to the 
FIFO Stack-Up Problem 

Since the FIFO Stack-Up problem is NP-complete in general, it is unlikely that there exist 
polynomial time algorithms for the problem. This motivates us to consider restricted versions of the 
problem, linear programming solutions, exponential time algorithms, parameterized algorithms, 
and approximation algorithms. 

5.1 Restricted Versions 

We have shown that the FIFO Stack-Up problem can be solved in polynomial time for the 
following special cases. 

• If the number k of sequences is constant then the FIFO Stack-Up problem can be solved 
in polynomial time 0{k ■ (A -|- 1)^), see Theorem 14.21 

• If the number p of stack-up places is constant we get a polynomial running time of 0{n + 
77 jP+2 ) by constructing the sequence graph, see Section [5.4.1l 

• If the number m of pallets is constant we get an algorithm with polynomial running time 

• ml), see Section [5.4.21 

• Finally, if the number n of bins is constant, then we restrict the input size, see Equation m 
and Remark 12.31 and therefore we get a constant running time 0{1). 

5.2 Linear Programming 

Next we give two linear programming approaches to solve the optimization version of the FIFO 
Stack-Up problem in which we have to compute a minimum number p such that a given list 
Q of sequences can be processed with at most p stack-up places. The first one computes a 
bin solution using 0{n'^) variables and the second one computes a pallet solution using 0{m^) 
variables. In Section |6] we will compare both solutions implemented in CPLEX and GLPK within 
an experimental study of running times. 

5.2.1 Computing a Bin Solution 

We have given a list Q oi k sequences and n bins 6 i,..., 6„ which use m pallet symbols. Our aim 
is to find a bin solution for Q, i.e. a bijection tt : [n] — [n] for the bins, such that by the removal of 
the bins from the sequences in the order of tt the number of needed stack-up places p is minimized. 

To realize a bijection tt we define binary variables xf G {0,1}, G [n], such that xf is 
equal to 1, if and only if bin bi is placed at position j by tt. In order to map every bin b^, i G [n], 
on exactly one position, i.e. to ensure tt to be surjective, we use the conditions 

n 

xi = 1 for every i G [n] 
i=i 
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and in order to map on every position j S [n] exactly one bin, i.e. to ensure tt to be injective, we 
use the conditions 

n 

xi = 1 for every j € [n]. 

Further we have to ensure that all variables xj, i,j S [n], are in {0,1}. We will denote the previous 
+ 2n conditions by PERMUTATiON(n, 

By TT the relative order of the bins of each sequence q G Q has to be preserved, consider for 
example the bin solution of the sequences of Example [221 The bins of the corresponding sequence 
are shown in black, the others are greyed out. 

ib5,be,bj,bs,bi,bg,b2,biQ,b3,b4) —>■ qi 
ib5,be,bT,bs,bi,bg,b2,bio,b3,b4) —>■ q2 

We have to consider the orderings of the bins given by the k sequences qi = ( 6 i,..., ),..., gj, = 

(fcnfc_i+i, • ■ ■, For every sequence qi, 1 < £ < k, and every bin bi of this sequence, i.e. 

rii-i + 1 < i < rii, we know that if bi is placed at position j, i.e. = 1 , then 

• every bin 6 ^/, i' > i of sequence qg, i.e. i < i' < ng, is not placed before bi, i.e. x\, = 0 for all 
j' < j, which is ensured by x^, < 1 — xf and 

• every bin bii, i' < i oi sequence qg, i.e. ng-i + 1 < F < i, is not placed after bi, i.e. x^, = 0 
for all j' > j, which is ensured by x^, < 1 — x^. 

Since we have 0{n'^) pairs {f ,j) and 0{N‘^) pairs (i', i), and since we have k sequences, there are at 
most 0{k ■ v ? ■ N"^) such conditions. We will denote these conditions by SequenceOrder(Q, xj). 
By an integer valued variable p we count the number of used stack-up places for some given 


sequence Q 

as follows. 

minimize p 

(5) 

subject to 


PERMUTATiON(n,x^), and SequenceOrder((3,x^) 

(6) 


and 

/(t, c) < p for every c G [tt- — 1] 

(7) 


and 

f{t,c)=( V ^ ( V ^i) 

(8) 


For the correctness note that subexpression f(t, c) is equal to one if and only if in the considered 
ordering of the bins there is a bin b' with plt{b') = t opened at a step < c and there is a bin b" 
with plt{b") = t opened at a step > c, if and only if pallet t is open after the c-th bin has been 
removed. 

Remark 5.1 Equation is propositional logic and not a linear function in standard form. 
Propositional logic can be reformulated for binary integer linear programming using the results 
ofm which show that every n-ary boolean function f{xi,... ,Xn) = Xn+i can be defined with a 
binary linear program using n-\-l boolean variables Xi,..., Xn+i ■ In order to express Equation m 
we need to define binary conjunctions and an n-ary disjunction. 

• Every conjunction xi A X 2 , xi,X 2 € {Oil} can be realized by introducing a new variable 
X 3 € { 0 , 1 } and three conditions 

X3 — xi < 0, X3 — X2 < 0, a:i -I- a;2 — Xa < 1 

such that finally X 3 = x\ A X 2 . 
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• Every n-ary disjunction V a ;2 V ... V a:„ can be realized by introducing a new variable 
Xn+i € { 0 , 1 } and n + 1 conditions 

Xi-Xn+1<0, X2-Xn+1<0, ■■■, Xn - Xn+1 < 0, Xi + ^2 + . . . + - X„+i > 0 

such that finally Xn+i = a;i V X2 V ... V 

Next these two transformations are applied on Equation m- We define m- (n — l) < binary 
variables g{t, c) G { 0 ,1}, t € {m\, c G [n — 1] such that 

9{t,c)= V xl, 

iG[n],j<c,plt{bi)—t 


which can be realized by 

for every i G [n],j < c,plt{bi) = t,t G [m],c G [n — 1] (9) 

and ( ^ x{) — g{t, c) > 0 for every t G [m],c G — 1] (10) 

Further we define m ■ (n — 1) < binary variables h{t,c) G {0,1}, t G [m], c G [n—l] such 
that 

Ht,c)= y xi, 

iG[n],j'>c,plt{bi)—t 

which can be realized by 

xl — h(t,c) < 0 for every i G [n],j > c,plt{bi) =t,t G [m],c G [n — 1] (11) 

and ( ^ x{)—h{t,c)>0 for every t G [m],c G [n — 1] (12) 

[n] ,j'>c,plt{bi)—t 

Finally we define rn • (n — 1) < binary variables f(t, c) G {0,1}, t G [m], c G [n — 1], such 
that fit, c) = git, c) A h{t, c), which can be realized by 

f{t, c) — gif, c) < 0 for every t G [m],c G [n — 1] (13) 

and f{t, c) — h{t, c) < 0 for every t G [m],c G [n — l] (14) 

and g{t, c) + h{t, c) — fit, c) < 1 for every t G [m] ,cG [n — l] (15) 

Theorem 5.2 For every list Q of sequences the integer linear program computes 

the minimum number of stack-up places p in a processing of Q. 

5.2.2 Computing a Pallet Solution 

By Theorem 13.21 the minimum number of stack-up places can be computed by the directed path- 
width of the sequence graph Gq plus one. In the following we use the fact that the directed 
pathwidth equals the directed vertex separation number |32) . 

For a digraph G = {V, A) on n vertices, we denote by n(G) the set of all bijections tt : [n] —>■ [n] 
of its vertex set. Given a bijection tt G n(G) we define for i G [n] the vertex sets L(i, tt, G) = {m G 
V I 7 r(n) < i} and R{i,TT,G) = {u G j 7 r(it) > i}. Every position i G [n] is called a cut. This 
allows us to define the directed vertex separation number for digraph G as follows. 

d-vsn(G) = min max |{u G L(i, tt, G) j G R{i, tt, G) : (v, u) G A}] 

7rGn(G) l<i<|V| 

An integer linear program for computing the directed vertex separation number for some given 
sequence graph Gq = (V, A) on m vertices is as follows. To realize tt we define binary variables 
x^ G {0,1}, i,j G [m], such that xf is equal to 1, if and only if pallet Vi is placed at position j by 
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TT. Additionally we use a variable w in order to count the vertices adjacent to the right side of the 
cuts. 


minimize w (16) 

subject to 

Permutation (m, Xj) (17) 

C 

and ^ Y (j, c) < w for every c G [m — 1] (18) 

1=1 

and Y{j,c)= \/ (xjAx^^,) (19) 

[m] ,Vi)€A 

For the correctness note that subexpression Y {j, c) is equal to one if and only if there exists 
an arc from a vertex on a position j' > c to a vertex on position j. 

Remark 5.3 Equation H19i) is propositional logic and not a linear function in standard form. In 
order to express Equation m we need to define binary conjunctions and an n-ary disjunction, 
see Remark \5.1[ 

We define m'^ binary variables X{i,i', j, j') G {0,1}; S [m], such that X{i,i',j,j') = 

x\ A xl,, which can be realized by 

Xih j) f) -xj <0 for every i, i!, j, j' G [m] (20) 

and X{i,i',j,j') — xj, <0 for every i,i', j, j' G [m] (21) 

and xl xl, — X{i,i', j, j') < 1 for every i,i',j,j' G [to] (22) 

We define 0{m^) binary variables Y{j,c) G {0,1} for j,c G [m — 1], j < c, such that 

Y{j,c)= V 

2,z^e[m] ,( 1 ;^/ ,Vi)CiA 

which can be realized by 

X{i,i,j,j') -Y{j,c) < 0 for every j' e {c+ 1,.. .,m},i,i' G [to], {vi',Vi) G A,j,cG [to - 1](23) 
and ( E (a i'J, /)) - y (j) c) > 0 for every j, c € [to - 1] (24) 

,Vi)CiA 

Theorem 5.4 For every list Q of sequences the integer linear program (0D-{Z2D,{ig)-® com¬ 
putes the minimum number of stack-up places p = w -\- 1 in a processing of Q. 

5.3 Exact Exponential Time Algorithms 

The running time of exponential time algorithms is often given in the (!l*-notation, which hides 
polynomial factors. 

Theorem 5.5 The FIFO Stack-Up problem can be solved in time 0*{2'^). 

Proof The directed pathwidth of a digraph G = (P, A) can be computed in time 0*12^^) by 
[5]. Since the sequence graph Gq can be constructed in time 0{n + k ■ m^) by algorithm Create 
Sequence Graph shown in Figured the statement follows by Theorem l3.2l □ 

□ 
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5.4 Parameterized Algorithms 

Within parameterized complexity we consider a two dimensional analysis of the computational 
complexity of a problem. Denoting the input by I, the two considered dimensions are the size |/| 
and a parameter k{I), see [9] for a survey. 

A parameterized problem is a pair (If, k), where 11 is a decision problem, I the set of all instances 
of n and K : I — >■ N is a so called parameterization or parameter. The parameter k{I) should be 
small for all inputs I £ I. Please note that the following running times are exponential, but for 
small values of k{I) they may be good in practice. 

• An algorithm A is an xp-algorithm with respect to k, if there are two computable functions 

f,g:N^N such that for every instance I G X the running time of A on / (with input size 
|/|) is at most /(k(/)) • A typical running time is A parameterized 

problem (If, k) belongs to the class XP and is called slicewise polynomial, if there is an 
xp-algorithm with respect to k which decides 11. 

• An algorithm A is an fpt-algorithm with respect to k, if there is a computable function 
/ : N —?► N such that for every instance I € I the running time of A on / (with input size 
|/|) is at most f{K{I)) ■ |/|‘^ for some fixed c G N. A typical running time is • |/p. 
A parameterized problem (11, k) belongs to the class FPT and is called fixed-parameter 
tractable, if there is an fpt-algorithm with respect to k which decides 11. 

Fxed-parameter algorithms have shown to be useful in several fields, among there are: phylo¬ 
genetics m, biopolymer sequences comparison [4], artihcial intelligence, constraint satisfaction, 
and database problems geometric problems [lOj . and cognitive modeling m- 

In order to show fixed-parameter intractability, it is useful to show the hardness with respect 
to one of the classes W[t] for some t > 1 which were introduced by Downey and Fellows [7] in 
terms of weighted satisfiability problems on classes of circuits. The following relations - the so 
called W-hierarchy - hold and all inclusions are assumed to be strict. 

FPT C W[l] C W[2] C ... C XP 

For the FIFO Stack-up problem we choose the number of sequences k or the number of 
pallets m and others as a parameter k(J) from the instance I, in order to obtain the following 
parameterized problem. 

Name: k(/)-FIFO Stack-up 

Instance: A list Q = {qi,..., qk) of sequences and a positive integer p. 

Parameter: k(/) 

Question: Is there a processing of Q, such that in each configuration during the processing of Q 
at most p pallets are open? 

We give two useful instruments to obtain fpt-results used in the next subsections. The first re¬ 
sult gives a connection between the existence of ILP formulations and fixed-parameter tractability 
w.r.t. the number of variables. 

Theorem 5.6 ([22],[25]) If for some problem 11 there is an ILP using £ variables, then 11 can 
be solved for every instance I in time 0{\I\ ■ 

The next remark states the running time of an exhaustive search for problems in NP. The 
correctness follows immediately by the definition of NP-completeness. 

Remark 5.7 ([31]) If some problem 11 belongs to NP, then there is some polynomial q such that 
n can be solved for every instance I in time 
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5.4.1 XP-Algorithms 

Parameterization by number of sequences k By Theorem 14.21 we obtain an xp-algorithm 
with respect to the parameter k. 

Theorem 5.8 There is an xp-algorithm that solves fc-FIFO Stack-Up in time 0{k ■ {N 1)^). 

Parameterization by number of stack-up places p The directed pathwidth of a digraph 
G = (F, A) can be computed in time 0(|A| • c by [55] and the 

sequence graph Gq can be computed in time 0{n + k ■ m?) which can be bounded by 0(n + m^) 
due to Remark 12.31 Thus the FIFO Stack-Up problem can be solved by Theorem 13.21 in time 
0{n + 

Theorem 5.9 There is an xp-algorithm that solves p-FIFO Stack-Up in time 0{n + 

5.4.2 FPT-Algorithms 

Parameterization by number of bins n We generate all possible bin orders (6.n-(i)! • ■ •: 
and verify, whether this can be a processing. This leads to a simple but very inefficient algorithm 
with running time 0{n^ ■ n\), where 0{n^) time is needed for each verification. 

Theorem 5.10 There is an fpt-algorithm that solves n-FIFO Stack-Up in time 0{n\ ■ n^). 

Alternatively we can apply Theorem 15.61 Remark 15.11 and Theorem 15.21 Integer program 
(ED-©,®-® has at most 4n^ + 1 G 0{n^) variables and a polynomial number of constraints. 
From Theorem 15.61 it follows that the FIFO Stack-Up problem is fixed-parameter tractable for 
the parameter n. 

Theorem 5.11 There is an fpt-algorithm that solves n-FIFO Stack-Up in time 0{\I\ ■ (4n^ -I- 

l)0{r^))_ 

Alternatively we can apply Remark 15.71 Since we assume that p < rn < n, the size of an 
instance I can be bounded by |/| G 0{n ■ log(n)), see Equation ©, and by Remark [5Tl it follows 
that the FIFO Stack-Up problem is fixed-parameter tractable for the parameter n. 

Theorem 5.12 There is an fpt-algorithm that solves n-FIFO Stack-Up in time 
for some polynomial q. 

Parameterization by number of pallets m Let (ti,... ,tm) be a permutation of the pallets 
of plts{Q). It can be checked in time 0{n ■ k) C 0{n'^) whether this permutation describes a 
processing of Q using at most p stack-up places, see Section 14.21 There are ml permutations of 
the pallets. If we enumerate all ml possible pallet orders and check each order, we can solve the 
FIFO Stack-Up problem. 

Theorem 5.13 There is an fpt-algorithm that solves m-FIFO Stack-Up in time 0(nf ■ ml). 

Alternatively we can apply Theorem l5.61 Remark l5.31 and Theorem l5.41 Integer linear program 
(fT6|-(p®.(f20|-([^ has at most -I- 2mf -I- 1 G 0{m'^) variables and a polynomial number of 
constraints. From Theorem 15.61 it follows that the FIFO Stack-Up problem is fixed-parameter 
tractable for the parameter m. 

Theorem 5.14 There is an fpt-algorithm that solves m-FIFO Stack-Up in time 0{\I\ ■ (m^ -\- 

2m^ + 

We also can apply the relation to directed pathwidth given in Theorem 13.21 Since the directed 
pathwidth of a digraph G = (F, A) can be computed in time 0(1.891'^!) by [53], and since the 
sequence graph Gq can be constructed in time 0{n -\- k ■ mf) by algorithm Create Sequence 
Graph shown in Figure jS] the next result follows by Theorem 13.21 
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Theorem 5.15 There is an fpt-algorithm that solves m-FIFO Stack-Up in time 0{n + k-w? + 

Alternatively we can apply the relation to directed pathwidth given in Theorem l3.2l and Propo¬ 
sition Let {Q,p) be an instance for the FIFO Stack-Up problem with values {n,m, k, N). 
We start to compute from Q the sequence graph Gq = {V, A) in time 0{n + k-rrL^). Then from Gq 
we build up the sequence system Qcq with values {n',m', k', N') in time 0(|F| + |A|) C 0{m?). 
By Theorem [32] and Proposition [331 list Q can be processed with at most p stack-up places, if and 
only if list Q' can be processed with at most p stack-up places. Further we know that k' € 0{m?) 
and N' = 2. By applying Theorem 14.21 for Q' we can solve the FIFO Stack-Up problem in time 
0{k' ■ {N' + 1)'=') C 0{m^ ■ 3™'). 

Theorem 5.16 There is an fpt-algorithm that solves m-FIFO Stack-Up in time 0{m‘^ ■ 3™^ -I- 
n + k ■ m^). 

Parameterization by combined parameters In the case of W[l]-hardness with respect to 
some parameter ^ a natural question is whether the problem remains hard for combined parameters, 
i.e. parameters (I', that consists of two or even more parts of the input. Since the existence of 
an fpt-algorithm w.r.t. parameter k is open up to now, we next conclude an fpt-algorithm with 
respect to parameter {k,m), by the result of Theorem 14.51 

Theorem 5.17 There is an fpt-algorithm that solves (fc,m)-FIFO Stack-Up in time 0{nf ■k'^). 

In practice k is much smaller than m, since there are much fewer buffer conveyors than pallets, 
thus this solution is better than 0{n^ ■ ml). In order to show a better fpt-algorithm with respect 
to parameter {k,m), we consider the result of Theorem 14.41 

Theorem 5.18 There is an fpt-algorithm that solves (fc,m)-FIFO Stack-Up in time 0{n'^ ■ {m-\- 
2 )'=). 

Since there are much fewer buffer conveyors than pallets in practice, this solution is better 
than 0{n^ ■ fc™). 

5.5 Approximation 

Since the directed pathwidth of a digraph G = (V, A) can be approximated by a factor of 
0(log^'® |F|) by [23], the FIFO Stack-Up problem can be approximated using the sequence 
graph Gq by Theorem 13.21 

Theorem 5.19 There is an approximation algorithm that solves the optimization version of the 
FIFO Stack-Up problem up to a faetor of m). 


6 Experimental Results 

Next we want to evaluate an implementation of algorithm Find Optimal Pallet Solution 
given in Figure [T3 and our two linear programming approaches given in (©-dZ]),®-® and (Hl])- 
(Ull),([201)-dill) using GLPK and CPLEX. 

6.1 Creating Instances 

Since there are no benchmark data sets for the FIFO Stack-Up problem we generated randomly 
instances by an algorithm, which allows to give the following parameters. 

• Pmax an upper bound on the number of stack-up places needed to process the generated 
sequences 
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• k number of sequences 

• m number of pallets 

• ^min and Tinax the minimum and maximum number of bins per pallet 

• d maximum number of sequences on which the bins of each pallet can be distributed 

The idea is to compute a bin solution B = ( 6 i,..., bn) with respect to the given parameters 
and to distribute the bins to the k sequences such that the relative order will be preserved, i.e. bi 
will be placed left of bj in some sequence, if i < j in B. The algorithm is shown in Figure [T5] 


Algorithm Random Instances 


#open := 0 ; open := 0 ; 

► 

CLVg := (rmin 4” I"max)/2 

► 

for all i := 1 to m step 2 do 

► 

r := random( 0 , rmax - ?’min) 

► 

no[t] := avg + r 

► 

no[i -1- 1 ] := avg — r 

► 

for all i G [m] and j £ [d] do 

► 

seq[i][j] := random(l,fc) 

► 


for all i e [k] do qi := {) 
unproc := [m]; i := 0 ; 
while i < n do 


Description of Functions and Variables: 
random(£, m): hit // choose some value in [£, m] 
at random 

no[l... m]: int // number bins for each pallet 
n: int // total number of bins, n = no[i\ 
seq[l... m\[l... d]: int // for each individual 
pallet there are up to d sequences, on 
which to distribute the bins 


i = i + 1 
if #open= p 

pit := choose some pallet of open at random 

else 

pit := choose some pallet of open U unproc at random 
if pit € unproc > first bin of chosen pallet 

^open := ^open+1; open := open U {pit}; unproc := unproc — {pit} 
r := random(l, d) 
s := seq[pZt][r] 

append bin bi to sequence qs 
no[plt] := no[plt]-l 

if no[plt] = 0 > last bin of pallet pit 

^open := ^open-1; open := open - {pit}; 


Figure 13: Construction of random instances for the FIFO Stack-Up problem. 


6.2 Implementations 

We consider the breadth first search based solution Find Optimal Pallet Solution for the 
FIFO Stack-Up problem given in Figure[T2l For several practical instance sizes the running time 
of this solution is too huge for computations, e.g. m = 254 and k = 10 lead a time of 0{n? ■ {m -I- 
2)^) = 0{'n? ■ 2®°). Therefore we used a cutting technique on the decision graph G = {V^A) by 
restricting to vertices v G V representing configurations Cq{v) such that ^open^Cq^v)) < Pmax 
and increasing the value of Pmax by 5 until a solution is found, see Figure fCT We have implemented 
algorithm Find Optimal Pallet Solution as a single-threaded program in C-l—|- on a standard 
Linux PC with 3.30 GHz CPU and 7.7 GiB RAM. 

Our two linear programming approaches given in (l5])-(r7)).(l^- (IT5|) and (fTOjl - lfTSll . ll^ - llMl) have 
been realized in GLPK v4.43 and GPLEX 12.6.0.0 and have been run on the same machine. GLPK 
is single-threaded, while CPLEX uses all 4 cores of the GPU. 
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Algorithm Iterative process 


Pmax ■— 5 5 

level := 0; > last level without cuts 

L[ 0 ] := ( ( 0 ,..., 0 ) ) O i[ 0 ]: level 0 only contains the initial configuration 

pred[ ( 0 ,..., 0 ) ] := 0 
val[ ( 0 ,..., 0 ) ] := 0 
open[ ( 0 ,... , 0 ) ] := 0 
while BFS(pmax) yields no solution do 
J^max ■— Pmax 5 


Figure 14: Iterative process of our cutting technique. 


Algorithm BFS(p) 


for £ := level +1 to m do > i: level in the decision graph 

for each configuration C G — 1] do \> L[P\: list of nodes in level I 

if {pred[C] ^ 0) 

Extend Path(C') 
for j := 1 to fc do 

O := open[C] U plt{bi.+i) 

for i := 1 to fc do > perform automatic transformation steps 

while the first bin in qi is destinated for an open pallet in O 

let Cs be the configuration obtained by removing the first bin of qi 
if {Cs is not in L[i] and ^open{Cs) < p) l> Cut 
append{Cs, L[i]) 
append{C, pred[Cs]) 
if {Cs is final) return Cs 

if £ — 1 7 ^ p then erase{C, L[£ — 1]) > remove C from list L[£ — 1] 

level: = p; 


Figure 15: BFS. 


6.3 Evaluation 

First we consider our implementation of algorithm Find Optimal Pallet Solution. In Table[2] 
we list our chosen parameters. For each assignment we randomly generated and solved 10 instances 
to compute the average time for solving an instance with the given parameters. Our results show 
that we can solve practical instances on several thousand bins in a few minutes. 

Next we consider our two linear programming models realized in GLPK and CPLEX. Since 
the size of the instances of Table [2] was so high that none of the ILP approaches was able to solve 
them, we generated much smaller parameters in the instances of Table [3l For each assignment we 
randomly generated and solved 10 instances in order to compute the average time for solving the 
same instances with the given parameters by our two linear programming models using GLPK 
and CPLEX. 

Our results show that algorithm Eind Optimal Pallet Solution can be used to solve 
practical instances on several thousand bins of the EIFO Stack-Up problem. Our two linear 
programming approaches can only be used to handle instances up to 100 bins and less than 10 
pallets. Since it is a commercial product of high licence cost CPLEX can solve the instances much 
faster than the open-scource solver GLPK and, as expected, the pallet solution approach is much 
better than the bin solution approach. 
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n 

PmsLX. 

Instance 

m k 

Pmin 

'?"inax 

d 

BFS and cutting 
Pallet 
Solution 

1 . 

1500 

14 

100 

8 

10 

20 

4 

0.14 

2 . 

1500 

14 

100 

8 

10 

20 

6 

0.12 

3. 

1500 

14 

100 

8 

10 

20 

8 

0.08 

4. 

2000 

14 

100 

8 

10 

30 

4 

0.17 

5. 

2000 

14 

100 

8 

10 

30 

6 

0.15 

6 . 

2000 

14 

100 

8 

10 

30 

8 

0.13 

7. 

2500 

14 

100 

8 

10 

40 

4 

0.22 

8 . 

2500 

14 

100 

8 

10 

40 

6 

0.10 

9. 

2500 

14 

100 

8 

10 

40 

8 

0.07 

10 . 

6000 

18 

300 

10 

15 

25 

5 

5.78 

11 . 

6000 

18 

300 

10 

15 

25 

7 

3.72 

12 . 

6000 

18 

300 

10 

15 

25 

10 

2.01 

13. 

7500 

18 

300 

10 

15 

35 

5 

7.78 

14. 

7500 

18 

300 

10 

15 

35 

7 

2.83 

15. 

7500 

18 

300 

10 

15 

35 

10 

2.07 

16. 

9000 

18 

300 

10 

15 

45 

5 

4.21 

17. 

9000 

18 

300 

10 

15 

45 

7 

2.23 

18. 

9000 

18 

300 

10 

15 

45 

10 

1.45 

19. 

12500 

22 

500 

12 

20 

30 

6 

92.52 

20 . 

12500 

22 

500 

12 

20 

30 

9 

52.74 

21 . 

12500 

22 

500 

12 

20 

30 

12 

42.81 

22 . 

15000 

22 

500 

12 

20 

40 

6 

103.24 

23. 

15000 

22 

500 

12 

20 

40 

9 

49.54 

24. 

15000 

22 

500 

12 

20 

40 

12 

32.85 

25. 

17500 

22 

500 

12 

20 

50 

6 

88.52 

26. 

17500 

22 

500 

12 

20 

50 

9 

38.61 

27. 

17500 

22 

500 

12 

20 

50 

12 

41.70 


Table 2: Running times in seconds for randomly generated instances for finding optimal solutions 
of the FIFO Stack-Up problem by algorithm Find Optimal Pallet Solution. 


7 Conclusions and Outlook 

In this paper we consider three graph models for the FIFO Stack-Up problem. Based on these 
models we have shown a breadth first search solution and two linear programming solutions to 
solve the problem. Further we have given parameterized algorithms w.r.t. several parameters 
which are summarized in Table 01 We also could give a first approximation result for minimizing 
the number of stack-up places. 

In our future work we want to determine the complexity of the FIFO Stack-Up problem for 
dg € {2,..., 5}. Further we intend to find better approximation algorithms, try to improve the 
running time of the given parameterized algorithms, and explore the existence of fpt-algorithms 
w.r.t. parameters k and p. Due Tamaki (Section 6 in |29| 1 the existence of an fpt-algorithm for 
the directed pathwidth problem w.r.t. the standard parameter is still open. By Theorem 13.21 such 
an algorithm would imply an fpt-algorithm for the FIFO Stack-Up problem w.r.t. p, and vice 
versa. 

We are also interested in on-line algorithms for instances where we only know the first c bins 
of every sequence instead of the complete sequences HIH]. Especially, we are interested in the 
answer to the following question; Is there a d-competitive on-line algorithm? Such an algorithm 
must compute a processing of some Q with at most p ■ d stack-up places, if Q can be processed 
with at most p stack-up places. First approaches for on-line algorithms for controlling palletizers 
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Instance 




GLPK 

GPLEX 









Bin 

Pallet 

Bin 

Pallet 


n 

Pmax 

m 

k r 

min 

^max 

d 

Solution 

Solution 

Solution 

Solution 

1 . 

15 

2 

3 

2 

4 

6 

2 

41.4 

0.1 

0.1 

0.1 

2 . 

20 

2 

4 

2 

4 

6 

2 

- 

0.1 

1.7 

0.1 

3. 

30 

4 

5 

4 

4 

8 

2 

- 

0.7 

66.8 

0.2 

4. 

48 

4 

6 

4 

6 

10 

2 

- 

2.5 

932.6 

1.2 

5. 

64 

4 

8 

5 

6 

10 

2 

- 

684.6 

- 

16.3 

6 . 

100 

5 

10 

5 

5 

15 

2 

- 

- 

- 

282.3 


Table 3: Running times in seconds for randomly generated instances for finding optimal solutions 
of the FIFO stack-up problem. Running times of more than 1800 seconds = 30 minutes are 
indicated by a bar (-). 


parameter 

k 

P 

m 

n 

{k, m) 

EPT 

7 

7 

+ 

-\- 

-k 

XP 

-f 

-f 

+ 

+ 

-k 


Table 4: Parameterized complexity of the FIFO Stack-Up problem 


are presented in [19]. 

In real life the bins arrive at the stack-up system on the main conveyor of a pick-to-belt 
orderpicking system. That means, the distribution of bins to the sequences, for example by some 
pre-placed cyclic storage conveyor, has to be computed. Up to now we consider the distribution 
as given. We intend to consider how to compute an optimal distribution of the bins from the 
main conveyor onto the sequences such that a minimum number of stack-up places is necessary to 
stack-up all bins from the sequences. 
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